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Abstract 

We study the neutral current flavour changing rare decay mode B —► X s + missing 
energy within the framework of theories with large extra spatial dimensions. The 
corresponding Standard Model signature is B —> X s + uu. But in theories with large 
extra dimensions, it is possible to have scalars and gravitons in the final state making 
it quite distinct from any other scenario where there are no gravitons and the scalars 
are far too heavier than the B-meson to be present as external particles. We give an 
estimate of the branching ratio for such processes for different values of the number of 
extra dimensions and scale of the effective theory. The predicted branching ratios can 
be comparable with the SM rate for a restrictive choice of the parameters. 
Keywords: Extra dimensions, Rare B decay 
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1 Introduction 

The investigations of flavour changing neutral current (FCNC) transitions of b-quark offer 
excellent opportunities to test the basic structure of the underlying theory. These processes 
are quite sensitive to QCD and possible long distance corrections as well as contributions 
from new particles in the loops. Consequently, such processes become useful tools for test¬ 
ing new physics as well. The measurement of B —> X s ~/ (and subsequently the exclusive 
channel B —> K*y) by CLEO |T] has been used to stringently constrain the parameter space 
of supersymmetric (SUSY) and various other theories (see for example |2j and references 
therein). The same is true for any other FCNC process. In particular, the quark level tran¬ 
sition b —> svu can turn out to be a very successful place for putting strong and meaningful 
bounds on the underlying theory. The Standard Model (SM) experimental signature in this 
case is the observation of the decay process B —> X s + missing energy in the form of neutri¬ 
nos. In any extension of the SM, decay of a b-quark to an s-quark and neutral, ultra light 
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particle(s) can mimic the SM process. 


It is well known that SM is plagued with the hierarchy problem and many possible solutions 
have been proposed in the past. However, the idea that the fundamental scale of gravita¬ 
tional interaction being quite distinct from the Planck scale ( rv./ 10 19 GeV) and possibly as 
low as O(TeV) has attracted a lot of attention. In the simplest version, as put forward by 
Arkani-Hamed et.al. (referred to as ADD model/scenario from now on) [3], the basic idea is 
the existence of n spatially compact large dimensions. The spacetime is a direct product of 
the four dimensional Minkowski space and the compact space spanned by the n spatial extra 
dimensions. The SM fields are all resticted to a 3-brane while gravity is free to propagate 
in the bulk. Seen from a four dimensional point of view, it simply means that apart from 
the usual massless graviton mode, there is a tower of Kaluza-Klein (KK) excitations from 
the gravity sector due to compactification. Also, if M* and Mpi denote the effective scale of 
gravity and the four dimensional Planck scale respectively and if R is the radius of compact¬ 
ification of the extra dimensions (assuming the same radius for all the n dimensions), then 
it turns out that they are related as follows 

M 2 Pl ~ R n M? +2 (1) 

Therefore, for large enough values of R , the effective scale M* can be as low as (P(TeV) and 
still be consistent with Eq. ©• Although the individual KK modes couple to the SM fields 
on the 3-brane by the ordinary gravitational strength but the presence of a large number of 
them makes the effective coupling appear to be of the order of TeV” 1 rather than Planck 
scale inverse. This means that gravity starts to become a strong force at TeV scales in sharp 
contrast to the situation in any ordinary theory of gravity. 


In the present study, we investigate the inclusive decay B —> X s + missing energy in the 
context of the ADD scenario. In the SM or most of its extensions, the missing energy is in 
the form of the neutrinos that escape without being detected and it is only new contributions 
to the relevant Wilson coefficients that are induced by interactions beyond SM. However, 
there is a sharp contrast in the case at hand. In the present scenario, there can be very light 
gravitons and associated scalars (the dilatons) that can be present in the final state. The 
situation is clearly distinct from SM or its usual extensions as in all such cases there are 
no spin-2 gravitons involved and the scalars, both neutral and charged, are far too massive 
compared to the b-quark to be present in the external legs. But this is possible in theories 
with large compact extra dimensions. It is precisely this advantageous aspect that we would 
like to exploit in the present study and see whether we get any meaningful results. 


2 b —> s + KK modes 

The effective Hamiltonian for b — > s transition in the SM is J3] 

n,„ = ~ v t y tb Y,c l (rio i (ri, 


( 2 ) 
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with 


0\ (^iCj)v-A (c-jbi) Y—A? 

O 2 — ( SiCi)y~A{Cjbj)v-Ai 

0 3 = (sibi) v _ A J2{Qj(lj)v-A, 
q 

0 4 = (^Iv-aE^v-a, 

q 

O5 {Sjbi)v—A ^ j ^Q.iQ.j)v+Ai ( 3 ) 

q 

Ofj = (s^v.aEMvaa, 
q 

°7 = J^Si a ^( m s p L + m b Ppt)biF^, 

and 

°s = ^Si(T^{m s P L + m b P R )T^G%. 

Apart from these there are semi-leptonic operators as well 


0(ee)y = (, Sibi) V - A (ee) v 


0(ee) A = ( Sibi) V -A(ee) A 


( 4 ) 


and finally 


0{yv) = (. Sibi) V - A (uv) V - A 


It is this last operator that is responsible for the relevant decay channel in SM. In order to 
calculate any process beyond these known operators, one will have to either write down all 
possible new operators respecting the symmetries of the low energy theory and then indi¬ 
rectly fix the associated coefficients or explicitly make the calculation. We follow the latter 
route. 


The coupling of the gravitons and the dilatons can be obtained from the low energy effective 
action UU- As expected, the gravitons couple to the energy-momentum tensor of the SM 
holds and the dilaton couples to the trace of the energy-momentum tensor. It is then straight 
forward to get the Feynman rules and compute the individual contributions. We consider 
the dilaton/scalar-KK emission first and comment on the analogous graviton emission later. 
The diagrams contributing to the dilaton emission process are shown in Fig.l. It may be 
important to mention here that: (i) there is no momentum dependent piece in the trace of 
energy-momentum tensor for the fermions and (ii) the dilaton-gauge boson-fermion-fermion 
vertex does not exist. However, both these are shown to be present in |5j. The fact that it 
is indeed so is not hard to see. The trace of the energy-momentum tensor for the fermionic 
plus gauge bosonic parts (including interaction term) is nothing but 

= ~ m fff + m v A2 ( 5 ) 
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In obtaining this, one needs to consider the interacting Heisenberg field equations rather than 
the free ones. The use of full interacting equations of motion for the fermion and the gauge 
boson results only in the mass terms for each of them (similar expression is obtained for the 
trace of the energy-momentum tensor in J7J) and thus there is no momentum dependence for 
the fermions involved in the Feynman rules. Also, the gauge boson-fermion- fermion-dilaton 
vertex will never be present (although a similar vertex for the gravitons will be there as in 
0). Thus, we differ on these points with the rules given in jjj and it may be important 
to emphasize again that in obtaining the energy-momentum tensor or its trace, nowhere 
has the on-shell condition been imposed. Therefore, the result Eq. © is exact and true in 
general. Similar remarks and results were quoted H1J in the context of a loop calculation in 
the Randall-Sundrum scenario, where the scalar in the theory (called the radion) couples, at 
the first order, with the trace of the energy-momentum tensor of the SM fields and is very 
similar to the dilaton in the ADD scenario as far as the couplings with the SM fields are 
concerned. There are less number of diagrams contributing to the dilaton emission process 
as compared to the graviton emission process. Also, in the limit of neglecting the strange 
quark mass, the diagram with dilaton being emitted from the s-quark line also vanishes. 


Clearly, all the diagrams are divergent and there is no underlying symmetry like the one 
present in SM that ensures cancellations between the individual diagrams to render a finite 
result. By naive power counting, one finds that there are hard divergences (see Appendix) 
which are quadratic, quartic and possibly even higher in the case of graviton emission. It 
now brings the question of handling such divergences in the context of effective field theo¬ 
ries. It has been very strongly advocated ID] that in any sensible effective field theory all the 
quadratic and higher divergences are cancelled by counter-terms arising out of the complete 
underlying high-energy theory. At the one loop level, it is thus the logarithmic dependence 
on the cut-off only that can be extracted from the low energy effective theory. Thus, it 
makes more sense to keep only the logarithmic pieces as far as the one loop calculation 
within an effective theory is concerned. We adhere to this approach in our calculations and 


thus use dimensional regularization throughout and at the end replace - by In 
the identification A = M*. 



with 


As mentioned earlier, it is expected that the presence of the scalars and gravitons in the 
external legs makes the situation more interesting as new operators are introduced. The 
operators are of the form (denoting the dilaton and the graviton fields by <F and h^ respec¬ 
tively) 

0(dilaton) ~ (s;(l ± 75 ) 64 )$ ( 6 ) 

and 

O(graviton) ~ (s* 7 M p bi ,(l ± 75 )b i )h 1 "' (7) 

(see Appendix for details) 

The invariant matrix element for the quark level process b —> s<I>^ for a dilaton of mass m n 
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is (the individual; contributions are listed in Appendix) 


M(b^ s$) 



( 8 ) 


where is the CKM factor, is the mass of the up-type quark in the loop. We have 
2 

neglected 0(-^f-) and higher terms in obtaining this expression. Using this matrix amplitude 

m w 

it is now trivial to compute the decay rate for the emission of a single dilaton of mass m n . 
For the inclusive process, it suffices to use the quark level amplitude to get the decay rate. 
One should sum over the tower of these dilatons till the b-quark mass scale. Following the 
summing techniques illustrated in |5], we get the following expression for the decay rate 




x 



Gl' 

2 

'1 

-b 

n 



1 

n + 4 




( 9 ) 


where we have made the replacement 



The expression for the graviton emission rate will look similar in its final form (see Appendix) 
but the rate itself is expected to be higher because of two reasons. Firstly, there are two 
extra diagrams (with the graviton being attached to either of the ferrmion-fermion-gauge 
boson vertex) and secondly, the graviton coupling does not have the u factor suppression 
present in the dilaton coupling. Thus, the rate is expected to be 

r (B -> X S G) ~ T T(B —> A. S T) (10) 

where T is a numerical factor expected to be of the order 10 or higher. However, it is 
instructive to go through the complete steps for the graviton emission process as well. The 
procedure is completely in analogy with the dilaton emission. The total rate is obtained by 
summing over various massive modes in the final state. Using the expression for the rate 
corresponding to a single graviton of mass m n , the summation essentially involves integrating 
over the tower of the massive modes till the b-mass scale. When this is done, it is found that 
there is a divergent behaviour for small n, the number of extra dimensions. This is like the 
soft emission rates diverging for some particular values of n and arises due to the structure 
of the integrand. We thus quote the values only for n > 4 (in particular n = 4, 6). In fact, 
the limit n —> 4 has also to be taken with care and in some sense in the spirit of analytic 
continuation or something similar. 
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3 Results 


The total contribution to the process D —> X s + KK modes is 

T(B -> X s $) + T(B X S G) 

Below we quote the branching fraction for the dilaton mode for different values of n and M*. 
The branching ratio for the graviton mode for n = 4 and n — 6 are as follows: 


Table 1: Br(B -> X, + $) 


n 

M* (TeV) 

Branching ratio 

2 

1 

1.28 x 10~ 3 

3 

1 

2.25 x 10“ 6 

4 

1 

4.93 x 10~ 9 

2 

5 

5.50 x 10~ b 

3 

5 

1.93 x 10~ 9 

4 

5 

8.45 x 10~ 13 

2 

10 

4.69 x 10~ 7 

3 

10 

8.23 x 10~ n 

4 

10 

1.80 x 10~ 14 

2 

50 

1.33 x 10~ 9 

3 

50 

4.68 x 10~ 14 

4 

50 

2.05 x 10~ 18 


Table 2: 5r(5 -»■ X s + G) 


n 

M* (TeV) 

Branching ratio 

4 

1 

1.02 x 10~ 4 

6 

1 

3.2025 x 10~ n 

4 

5 

1.72 x 10~ 8 

6 

5 

2.2 x 10~ 16 

4 

10 

3.67 x 10~ 1U 

6 

10 

1.17 x 10~ 18 

4 

50 

4.17 x 10~ 14 

6 

50 

5.33 x 10 -24 


These numbers are to be compared with the SM expectation, Br(b —> svv) = 5 x 1CB 5 
and the experimental limits on the same, Br(b —> svv) exp < 6.2 x 1CT 4 at 90% confidence 
level PI- From the Table 1 it is clear that for n = 2 and M* = 1 TeV, the contribution 
far exceeds the SM prediction and the experimental limits and should have been observed. 
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Moreover, one encounters some divergent results for the graviton emission process for smaller 
values of n which may be looked upon as arising due to something similar to gravi-strahlung 
typical to these smaller values of n. We therefore cannot say anything very satisfactorily 
for these values unless a more careful and detailed analysis is done. Thus we restrict our¬ 
selves to dilaton rates for these values of n and don’t consider graviton rates for those values. 


Also, evident is the fact that for higher values of n and/or M*, the branching fraction is 
smaller. The constraints obtained from the collider experiments m allow even smaller val¬ 
ues of M* than quoted, thus allowing for the possibility of some more enhancement in the 
branching ratio predictions. However, strong constraints from the supernova SN1987A ob¬ 
servations m and cosmology na rule out such choices of n and M*. For n = 2, SN1987A 
constraints imply M* ~ 50 TeV. Even if one takes care of all the uncertainties entering the 
supernova calculations, M* = 1 TeV seems to be far from reachable. The other values quoted 
in the table are more or less allowed by these constraints but seem to be contributing very lit¬ 
tle. However, for n — 4 and M* greater than a few TeV, the sum of the dilaton and graviton 
emission rates can be sizeable and may compete with the SM rate. The constraints coming 
from other astrophysical and cosmological processes dmu like cosmic diffuse gamma ray 
background, early matter domination etc are even stronger than quoted above and rule out 
any chances of extra dimensional contribution to the decay process being comparable to 
the SM rate and thus being observable at the future B-factories. Also, given such strong 
constraints, the idea of having different radii for different compact dimensions, or atleast 
partially, does not seem to be of much significance. 


It may be useful to mention here that the mode B —> X S G can contribute to SM process 
B —* X s vv via G —> vu and the sum of the SM and the extra contribution will form the 
complete matrix element. But, in the case of the scalar mode, B —» V s <h, there can be no 
cascade decay of <h to vv. Hence, until and unless the produced <h decays into low mass 
particles and escapes as missing energy, the energy and angular spectrum of the hadronic 
junk, X s , will be markedly different from the corresponding distribution observed in the case 
of B —> X s vv coming from any theory. The same is true for the graviton emission process 
provided it does not decay further into vv or any other pair of low mass particles. From the 
structure of O(graviton) it is evident that the operator corresponding to b —> svv will have 
the form 

O(vv)graviton ~ (sibi) v ± A (vv)v 

which is pure vector current for the vu pair, a feature different from the SM again where the 
intermediate particle responsible is the Z-boson, making the angular distribution different. 
Also observe that if instead of following the approach advocated in jUJ, we had retained 
quadratic and higher divergences, we would have obtained enormously large, intolerable 
numbers. This again justifies the retaining of only logarithmic terms and in turn justifies 
our use of dimensional regularization in carrying out the calculations, which otherwise would 
have to be carried out with a hard ultra-violet cut-off. 


In conclusion, we can say that for two extra dimensions and with M* > few TeV (as 
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marginally allowed by the constraints from astrophysics and cosmology), the desired branch¬ 
ing ratio can be sizeable. Other choices of n seem to be too small to compete with the SM 
rate. Also, for n > 2, the radius of compactification, R , becomes too small for its effects to be 
probed at future mechanical experiments testing the gravitational law at small length scales 
(eg. see [T3] for details of these experiments). Thus, theories with large extra dimensions 
predict a competing value of the branching ratio for the process B —» A s + missing energy 
(in the form of soft dilatons and gravitons) for a very restrictive choice of the parameters. 
One can thus hope that future observation of the decay mode B —> A" s + missing energy, 
along with the corresponding exclusive modes, will be able to testify theories with large extra 
dimensions and if not completely ruled out, yield severe constraints on the number of extra 
dimensions and the effective scale of gravitational interactions in these theories. 


Appendix 


We outline the steps involved in the evaluation of a typical loop diagram and quote the 
expressions for individual matrix elements. For the sake of illustration of the calculation, we 
evaluate the loop diagram (b), where the dilaton is emitted from the W-propagator. We call 
this matrix element M.^ b y Label the momenta as p bl p s and corresponding to b-quark, 
s-quark and the dilaton and denote the loop momentum variable by k. The matrix element 
can be written as 


M (b) = 


dk- 


[v 


fia 




][rfP 


k v kP 1 

<m I 


X 


(k 2 - m%y)[(k - p$) 2 - m 2 w ][(p b - k) 

Vap s(p*)[7mu(l - 75) (A- h + - 75 )]6(p&) 


m- 


( 11 ) 


where 




2 

wmi ; 


dk = 


d A k 

dt ) 1 


In the above expression, g is the weak coupling constant, G; is the relevant CKM factor, 
is the mass of the up-type quark in the loop, k 2 = 16ttGn gives the gravitational coupling 
in terms of Newton’s constant GV and w is a number that depends only on the number of 
extra dimensions, n, and is given by 


u 


2 


2 

3(2+ n) 


Clearly, u < 1 and gives a supression factor. Also, it can be seen very clearly from the 
above expression for the matrix element that the ^-integral is not logarithmically divergent 
by power counting but there are higher divergences present. The reason for the presence of 
such hard divergences is the fact that there is no underlying symmetry, like the one present 
in the SM, to ensure nice cancellations and render a finite or at the most logarithmically 
divergent result. This can be attributed to the non-renormalizable nature of these interac¬ 
tions. The calculation has been done in tell unitary gauge and it may seem that these power 
divergences are arising due to the form of the massive gauge boson propagator. However, 














it should be noted that if instead one works in some other convenient gauge, say Feynman 
gauge, the massive gauge boson propagators have a more well behaved high energy structure 
but there are additional diagrams involving the asociated ghosts which essentially are the 
badly behaved longitudinal components of these massive gauge bosons and also that their 
interaction vertices contain momentum dependence. Therefore, with some effort, various 
terms in the calculation in some other gauge can be combined to give the same results. This 
is also evident from the Feynman rules in [5] where the vertices have an additional gauge 
dependence. But we follow the thumb rule || that only logarithmically divergent terms 
be retained and thus use dimensional regularization (in the full spirit of an effective theory 
calculation) throughout our calculation. The integral can be evaluated using standard proce¬ 
dure of Feynman parameterization, shifting the variables and carrying out the integrations. 
We don’t evaluate the finite pieces and only retain divergent terms. After the Feynman 
parameterization and shifting of momentum variables, the terms that are not multiplied by 
rrii (or its powers) drop out due to the uuitarity of the CKM matrix (the GIM mechanism). 

Retaining just the various divergent terms in the integral and employing dimensional reg- 

2 

ularization (neglecting terms £?(-%-) and higher), we arrive at the following expression for 

m w 

the matrix element corresponding to the diagram (b): 


M (b) = 


y 167T 2 
s(Ps) 




m? 


2V2, 


UJK- 


m 


w 


( 12 ) 


3 2 

-m s ( I- 75 ) + -m b ( I + 75 ) 


KPb) 


where rn b and m s are the masses of the b- and s-quark respectively. Also, using the relation 
v§ = gTT" we can ex P ress the result in terms of Fermi-constant Gp. 


The expressions for other diagrams are as follows: 


M (a) = 


167T 2 


4£& 


9 


,2y/2, 


COK 


™i 

rri'w 


3 


s(Ps) 


m s ( I- 75 ) + m b ( I + 75 ) 


KPb 


(13) 


Al( c ) = 


167T 2 




0 

mf 


UJK- 


m 


s(p B ) 


2^2 

m b (l - 75 ) - m s (l + 75 ) 


w 


3 


m b m s 


e (mi — m- 


KPb ) 


(14) 


M-(d) 




2 

iOK, 


jn 



(1 \ m b m s 

V e / ( m l - m i) 


s(Ps) 


m b (l - 75 ) + m s (l + 75 ) 


KPb 


(15) 
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Adding all these contributions we arrive at Eq. ( 1211 - 


Similar to the dilaton emission process, we take a look at the graviton emission process and 
the structure that is obtained. We compute the graviton emission process corresponding to 
diagram (a). Following the same steps as outlined above we get for the matrix element and 
neglecting rn s 


M 


qraviton 

(a) 



s(Ps)laPb^ ~ 75 )KPb) C 


(16) 


where e a/3 is the polarization tensor for the graviton field. Other diagrams give similar 
contributions. Due to momentum dependent factors in the Feynman rules (when a graviton 
is involved), the divergence structure is even worse and naive power counting shows that 
divergences higher than quadratic are present. We don’t work in the cut-off scheme and 
the use of dimensional regularization again ensures that we retain only the logarithmic 
divergences. Note that there is no factor of uj appearing in the above expression and thus 
the supression due to this factor is absent in the graviton emission process and depending 
on the number of extra dimensions etc., the rate for the graviton emission process can be 
larger than the dilaton emission rate by a factor of 10 or more. This extra multiplicative 
factor has been identified as T in Eq. m- The graviton emission rate (for a single graviton 
of mass m n ) is 


T(B -> X s G {n) ) = 


X 



In 




m 2 n f{2ml + 3 m 2 n ) 

"4 



(17) 


Summing over all the states lying till the b-quark mass scale gives the final expression for 
graviton emission rate. Clearly, the form for the expression is similar to that of the dilaton 
case and there are no factors of u present. It is important to remark that from the above 
expressions, it is clearly evident that the oerator structure for the dilaton or graviton emis¬ 
sion processes is indeed the ones quoted in Eqs. ©-©• 
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Figure 1: The diagrams contributing to the dilaton emission. We label the diagrams as (a), 

(b), (c) and (d) moving clock-wise from top left. For the graviton emission there are two 

more diagrams with the graviton being hooked to either of the two ends of the loop. 
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